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ABSTRACT 


Since  the  initiation  of  electron  radiation  therapy, 
most  of  the  experimental  efforts  to  understand  the 
transport  of  electrons  in  a  medium  have  been  concentrated 
on  the  investigation  of  dose  deposition  by  electrons. 
Because  the  dose  distribution  of  an  electron  beam  is  only 
the  end  result  of  the  transport  of  electrons ,  it  has  not 
been  possible  to  predict  the  dose  distribution  in  a  complex 
medium.  In  addition,  the  prediction  of  a  dose  distribution 
even  for  a  simple  medium,  such  as,  an  all  water  "phantom", 
cannot  be  made  generally  because  of  the  the  degradation  of 
the  accelerated  beam  caused  by  the  method  of  shaping  a 
broad  beam(e.g.,  scattering  foil  and  collimator). 

In  this  study  of  electron  transport  in  a  homogeneous 
semi-infinite  medium,  the  linearized  Boltzmann  transport 
equation  is  used  with  the  boundary  condition  of  an 
individual  pencil  beam  is  used.  To  facilitate  the 
mathematical  treatment  of  the  problem,  the  medium  is 
divided  into  two  regions.  The  first  region  is  where  the 
path-length  of  an  electron  is  still  far  shorter  than  the 
"diffusion  length"  of  the  medium.  By  applying  the 
"Fokker-Plank"  approximation  we  obtain  the  angular 
distribution  of  electrons  and  the  average  lateral,  and 
longitudinal  displacements.  The  second  region  is  where  the 
orientations  of  electrons  has  become  very  diffuse.  In 
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order  to  take  care  of  the  energy  dependence  of  the 
scattering  cross-section,  we  use  the  "continuous  slowing 
down  approximation".  The  angular  density  function  is 
expanded  in  terms  of  spherical  harmonics  and  the  scattering 
cross-section  is  also  expanded  by  the  Legendre  polynomial. 
The  series  is  truncated  at  m=l  and  n=l.  The  spatial 
distribution  of  electrons  is  obtained  by  solving  the 
resulting  coupled  equations. 


Results  are  compared  to  those  obtained  by  Fermi  and 
discrepancies  are  found  in  angular  distributions  of 
electrons.  Furthermore,  it  is  pointed  out  that  Fermi’s 
solution  has  a  serious  drawback  because  the  depth  of  the 
medium  is  used  as  an  energy  indicator  instead  of  the  actual 
path-length  of  eletrons.  Its  application  to  the  radiation 
transport  problem  is  thus  inappropriate,  despite  its 
popularity  in  applications  to  medical  physics  problems. 
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CHAPTER  I 
INTRODUCTION 


1.1  History  of  Radiation  Therapy  with  Electrons 

In  1934,  Brasch  and  Lange(l)  reported  on  the  use  of 
"fast  cathode  rays"  as  a  method  of  radiation  therapy.  This 
is  apparently  the  first  paper  dealing  with  the  use  of 
external  beams  of  electrons  as  a  therapeutic  agent.  They 
cited  as  advantages; 

(a)  well  defined  and  adjustable  range  of  penetration; 

(b)  ability  to  deliver  more  dose  at  depth  relative  to  the 
dose  at  surface  ("skin  sparing"); 

(c)  enhanced  biological  response;  and 

(d)  the  possibility  of  magnetic  directabili ty . 

The  energy  of  the  electrons  they  used  was  between  0.7  MeV 
and  2  MeV,  with  a  range  of  up  to  approximately  1  cm  in  soft 
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tissue.  The  first  mention  of  the  use  of  higher  energy 
electron  beams  was  made  by  Hass  et  al.(2)  in  1954.  These 
authors  pointed  out  the  advantages  described  previously  and 
also  stated  that  the  skin  reaction  was  not  worse  than  that 
experienced  with  other  types  of  radiation.  This  paper 
dealt  with  electron  beams  produced  with  a  22  MeV  betatron. 

Electron  radiation  therapy  was  initiated  mainly 
because  of  the  advantageous  penetration  characteristics (a 
and  b)  and  it  evolved  rapidly  with  the  development  of 
reliable  higher  energy  accelerators.  The  early  method  of 
dose  administration  was  based  largely  on  the  energy 
deposition  data  measured  in  a  homogeneous  medium,  such  as  a 
water  "phantom" ( Figure  1).  This  approach  is  based  on  the 
assumption  that  the  volume  of  the  patient's  body  to  be 
irradiated  is  of  unit  density  and  water-like  in  atomic 
composition.  The  dose  distribution  of  an  electron  beam  is 
assumed  not  to  be  severely  distorted  as  it  penetrates 
heterogeneous  tissue.  This  dose  administration  method  is 
therefore  based  on  the  simple  consideration  of  the 
cross-sectional  area  and  the  depth  of  the  tumor  from  the 
skin  surface,  which  determine  the  field  size  and  beam 
energy,  respectively. 

In  1961,  Ovadia  and  McAllister ( 3 )  reported  improved 
"skin  sparing"  and  the  reduction  of  dose  to  normal  tissue 
overlying  tumors  with  the  use  of  a  "grid"  placed  in  the 
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Fig.  1.  Depth  dose  curves  of  broad  beam  of  electrons 
of  various  energies  incident  on  a  water 
phantom 
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incident  beam.  The  grid  was  constructed  of  a  lead 
mask(1.25  cm  thick)  with  holes  of  various  size.  These 
authors  optimized  the  hole  size  and  spacing  to  achieve  dose 
reduction  in  the  skin  and  "transit"  tissue.  In  1963,  there 
was  another  report  of  a  similar  effect  without  the  use  of  a 
grid  by  Carpender  et  al.(4)  who  described  electron 
radiation  therapy  with  a  scanned  "pencil  beam".  This 
concept  of  a  "pencil  beam"  is  important  and  will  be 
emphasized  in  this  thesis.  They  found  less  skin  reactions, 
and  gave  several  explanations ; 

(a)  "grid-effect"  caused  by  the  raster  scan  pattern, 

(b)  higher  local  dose  rate, 

(c)  uncontaminated  beam  of  electrons  achieved  by  not  using 
a  scattering  foil  to  produce  a  macroscopic  beam  by 
spreading  a  pencil  beam. 

However,  they  were  doubtful  about  the  first  explanation, 
because  the  scanning  was  performed  continuously  to  give  a 
homogeneous  surface  dose.  Furthermore,  daily  variations  in 
positioning  of  the  field  on  the  patient  would  further 
diffuse  or  blurr  the  raster  pattern. 

Another  significant  finding  was  the  actual 
alteration  of  the  depth  dose  curve(Figure  1),  caused  by 
tissue  inhomogeneity.  Taply  and  Fletcher(5)  summarized 
various  clinical  complications  caused  by  overdosage  due  to 
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the  presence  of  lung  in  1964.  Several  years  before  this, 
simple  corrections ( i . e . ,  CET)  for  tissue  inhomogeneity  had 
been  developed  by  Laughlin.(6)  Basically,  this  method 
relies  upon  the  modification  of  standard  central  axis  dose 
data  on  the  basis  of  the  average  density  of  tissue  actually 
traversed.  At  that  time,  the  lack  of  information  of  the 
anatomy  of  the  patient  presented  an  unsolved  problem  and 
more  complex  schemes  which  consider  electron  transport  were 
not  developed.  With  the  advent  of  X-ray  computed 
tomography  which  provides  not  only  the  geometry  of  the 
patient’s  anatomy  but  also  an  iri  vivo  map  of  tissue 
densities,  one  of  the  fundamental  obstacles  to  accurate 
treatment  planning  was  removed .( 9 , 11 , 12)  However,  in  spite 
of  this  remarkable  progress  in  obtaining  patient-specific 
data,  the  complexity  of  electron  transport  in  a 
heterogeneous  medium  has  yet  to  be  fully  understood  and 
this  remains  as  a  problem  today. 


1.2  Physical  Description  of  Broad  Electron  Beams 

1.2.1  Beam  Energy 

The  energy  is  perhaps  the  most  important  parameter 
in  that  it  is  directly  related  to  the  location  of  abrupt 
fall-off  in  the  depth  dose  curve,  as  measured  in  a 
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homogeneous  medium.  Selection  of  an  appropriate  energy 
permits  the  treatment  of  a  lesion  located  near  a  surface 
and  extending  to  a  specified  depth  (e.g.  chest  wall).  In 
fact  this  trend  of  depth  dose  curves  to  fall  off  sharply 
has  been  used  reciprocally  to  define  the  "mean  energy"  of  a 
broad  beam  of  electrons. 


The  nominal  energy  of  electrons  is  fixed  and 
uniquely  defined  by  the  accelerating  potential.  Even 
though  there  is  a  small  spread  in  energy(<l%)  within  the 
accelerator,  we  can  characterize  this  energy  by  a 
monoenerget ic  value  E  ,  the  "accelerator  energy".  The 
accelerated  narrow  beam  of  electrons,  however,  has  to  pass 
through  exit  windows  of  the  accelerator  space  ,  scattering 
foils,  mirrors,  dose  monitor  chambers,  and  air  layers  of 
considerable  total  thickness  before  reaching  the  body 
surf ace . (Figure  2)  A  cascade  of  energy  losses  in  these 
materials  shifts  the  electron  energy  to  less  than  E  , 

3. 

and  also  broaden  the  spectrum  as  a  result  of  fluctuations 
in  energy  loss.  Provided  these  fluctuations  in  energy  loss 
are  very  small,  one  can  still  speak  of  a  single  energy 
value,  Eq  ,  at  the  body  surface; 


AE 


tot 


where  Ae  t  total  mean  energy  loss  in  the  layers 

traversed  by  the  beam  before  reaching  the  patient.  In  a 
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ELECTRON  BEAM 


I 

i 


Fig.  2.  Typical  treatment  head  of  medical 

accelerator  used  for  electron  radiation 
therapy 
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practical  situation,  this  quantity  is  difficult  to 
calculate,  although  a  measurement  is  possible ( 49 ) .  Beyond 
layers  of  the  tissues  of  appreciable  thickness  in  an 
irradiated  body,  one  has  to  consider  a  further  distribution 
of  the  electron  fluence  in  energy,  which  is  influenced  by 
the  incident  spectrum,  further  energy  losses  in  penetrating 
into  the  body,  and  the  generation  of  secondary  electrons. 
The  distribution  may  be  characterized  by  the  most  probable 
energy,  the  mean  energy,  and  the  maximum  energy  at  a 
particular  site  in  the  body. 

According  to  the  ICRU(17),  there  are  four  possible 
statements  of  energy  currently  in  use; 


(a)  E  ,  with  description  of  the  scattering  foils(if 

cl 


any)  ; 

(b) 

the  "most 

probable  energy"  at  the  surface  of  the  body  ; 

(c) 

the  "mean 

energy"  at 

the  surface  of  the 

body  ; 

(d) 

energies 

determined 

experimentally 

using  various 

methods 

that  yield 

values  of  one  of 

the  above  three 

energies . 

These  quantities  are  shown  in 

Figure  3. 

$E( arbitrary  units) 
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Phantom 


\  Window  Accelerator 


Fig.  3.  Definition  of  various  energies  used  in  electron 
radiation  therapy.  $  is  the  distribution  of 
fluence  in  energy.  (17) 


E^:  The  accelerator  energy 

E  :  The  maximum  energy  at  the  surface 
max 

Eq :  The  energy  at  the  surface  of  medium 

E  :  The  energy  at  depth  Z 
z 
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1.2.2  Energy  Spectrum 

The  quality  of  the  electron  beam  is  characterized 
fully  by  its  energy  spectrum,  jr  (assuming  no  photon 
contaminations).  The  energy  spectrum  determines  the  yields 
of  energy-dependent  physical ,  chemical ,  and  biological 
effects.  The  energy  spectrum  also  has  a  considerable 
clinical  significance  because  it  influences  the  actual  dose 
distribution  of  electron  beams  (e.g.  sharpness  of  the  dose 
fall-off  region).  Moreover,  the  stopping  power  of  the 
medium,  which  is  energy-dependent,  is  the  basis  of  the 
calculation  of  the  absorbed  dose.  In  defining  the 
spectrum,  the  interesting  quantity  is  the  particle  fluence, 
that  is,  the  number  of  electrons  passing  through  a  unit 
area,  per  unit  time,  per  unit  energy(17).  Experimental 
information  about  the  fluence  within  an  irradiated  medium, 
is  very  hard  to  obtain,  often  one  must  rely  on 
calculation. (18, 19) 

Another  important  quantity  related  to  energy  spectra 
of  electrons  is  the  distribution  of  electrons  per  unit 
solid  angle  for  a  particular  direction  of  motion.  This 
quantity  is  essential  information  for  understanding 
influence  of  tissue  inhomogeneity  on  electron  dose 
distributions.  Multiple  scattering  of  electrons  in  a 
complex  medium  poses  tremendous  difficulties  in  calculating 
or  measuring  this  fundamental  quantity. 
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1.2.3  Field  Size 

Because  of  clinical  requirements,  the  field  size  of 
the  external  electron  beam  must  be  variable.  As  part  of 
the  specification  of  the  electron  beam,  the  field  has  to 
have  a  shape  and  size  compatible  with  the  spatial 
distribution  of  the  tumor  and  surrounding  healthy  tissue. 
Generally,  rectangular  fields  are  produced  by  the 
collimators  in  the  treatment  head(Figure  2),  but  these  can 
be  shaped  further  with  secondary  block  to  produce 
irregularly-shaped  fields. 


1.2.4  Field  Uniformity  And  Symmetry 

Field  uniformity  and  symmetry  are  specified  at  a 
fixed  depth  in  a  homogeneous  phantom  for  a  portion  of  the 
geometric  field  area  at  that  depth.  The  reference  depth  is 
often  taken  as  the  depth  along  the  central  axis  at  which 
90%  of  the  maximum  dose  is  delivered.  However,  if  we  are 
truly  interested  in  understanding  field  uniformity  and 
symmetry,  we  should  also  consider  the  electron  fluence(the 
number  of  electrons  per  unit  solid  angle  for  a  specific 
direction  of  motion,  per  unit  energy,  per  unit  volume). 
The  absorbed  dose  is  only  the  end  product  of  interactions 
between  the  electrons  and  the  constituent  atoms  of  the 
medium.  A  uniform  absorbed  dose  distribution  over  the 
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field  can  be  formed  by  many  different  fluence 
distributions,  which  are  more  of  fundamental  importance. 


1.2.5  Mass  Stopping  Power 

An  understanding  of  the  interaction  of  charged 
particles  with  matter  is  important,  because  the  detection 
and  influence  of  all  types  of  radiation  on  all  subtances, 
living  or  not,  is  through  this  interaction.  For  example, 
the  detection  of  neutral  particles (e  .  g .  neutrons,  photons) 
can  only  take  place  indirectly,  i.e.,  after  these  have 
produced  secondary  charged  particles (e . g . ,  electron, 
proton)  which  dissipate  enough  energy  to  be  detectable. 
Because  several  processes  contribute  to  the  slowing  down  of 
a  fast  charged  particle,  it  is  convenient  to  discuss 
stopping  power  in  three  different  energy  intervals; 

-4  2 

(a)  low  energy,  below  approximately  10  ire  (  m  is  mass  of 
electron ) 

(b)  intermediate  energy 

2 

(c)  high  energy,  above  nc  for  all  particles. 

For  the  low  energy(a),  elastic  collisions  between 
the  particles  and  "whole"  atoms  of  the  medium  are 
important.  In  the  intermediate  energy  region,  the 
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principal  interaction  is  between  the  particle  and  the 

atomic  electrons  of  the  stopping  medium.  Electrons  are 

excited  into  higher  quantized  energy  states ( exci tat ion )  or 

continuous  energy  states ( ionization ) .  This  inelastic 

collision  process  is  well  understood  through  Coulomb  forces 

and  the  stopping  power.  The  average  energy  loss  per  unit 

length  of  path  can  be  computed  from  a  knowledge  of  a  few 

parameters  measured  experimentally.  The  bulk  properties  of 

the  material  (e.g.  density)  also  influence  the  interaction 

process  and  tend  to  decrease  the  stopping  power.  This 

phenomenon  is  called  the  "density"  effect  or  "polarization" 

effect.  Nuclear  interactions  occur  only  rarely(  1 
g 

per  10  electron  interactions),  but  their  large  effect 
(i.e,  severe  energy  loss)  is  of  importance  for  the  upper 
range  of  intermediate  energies(b)  of  heavy  charged 
particles . 

For  high  energies(c),  radiative  interactions  or 
"bremsstrahlung"  is  of  major  importance. 

The  interactions  of  electrons  have  two  distinctive 
features;  (a)  For  each  collision  with  atomic  electrons, 
large  energy  losses  can  occur;  and  (b)  electrons  with 
energies  of  even  only  a  few  hundred  keV  will  be  influenced 
by  relativistic  effects.  In  general,  it  is  assumed  that 
after  an  electron-electron  Coulomb  collision  the  emergent 
electron  with  the  greater  energy  will  be  considered  the 


14 


"primary" .  The 

^  2 

o  f  MeV/cm  g 

stopping  power; 


mean  energy  loss,  expressed  in  units 
,  is  given  by  the  collisional  mass 


Q.)5ig\ 

P1  ) 


Z 

A 


2(T  +2) 

(i/mc2-)1. 


■+*  F  [ly&)  ~  S' 


where  4”  2  T  and 

F  =  -\-0‘+Jln[(T-A}A]  +  [l/(T-A)] 

4-A*>  (2.T+ 1 ) 

4-  - - — - - - — - — 

(T+I)a 


Here  "X  =:  (T/mc2)  >  d  is  the  density  (or  polarization) 

correction,  and  me  =511,006eV;  A  is  the  maximum  energy 

2 

transfer  of  interest  in  units  of  me  .  The  above 
equations  were  written  containing  A  explicitly  so  that 
they  would  be  useful  in  the  description  of  the  linear 
energy  transfer (LET )  or  restricted  stopping  power (45). 

A  second  important  effect  for  the  energy  loss  of 
electrons  is  the  emission  of  photons  due  to  acceleration 
(or  decelleration )  of  a  moving  charge.  This  bremsstrahlung 
energy  loss  increases  rapidly  with  an  increase  in  energy  of 
the  electrons  or  atomic  number  of  the  absorber.  A  critical 
energy,  ,  can  be  defined  at  which  the  energy  loss  due  to 
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Coulomb  collision  is  equal  to  loss  due  to  bremsstrahlung ; 

Tc  ^  [700/ (Z  +  1.2)]  MeV. 

For  human  soft  tissue,  Tc  is  approximately  85  MeV.  An 
estimate  of  the  ratio,  r,  of  the  energy  loss  due  to 
bremsstrahlung  to  that  due  to  collisions  is  given  by 

r  ^  (  TZ/700  )  . 

Thus  for  soft  tissue(Z^7)  and  electron  beams  up  to  30  Mev , 
r  ^>30%.  For  "hard"  tissue ( Z ^ 14 ) ,  such  as  bone,  this 
proportion  rises  to  ^ 60  %. 


1.2.6  Mass  Angular  Scattering  Power 

When  a  charged  particle  traverses  a  slab  of  finite 
thickness,  it  undergoes  a  large  number  of  collisions,  most 
of  which  produce  very  small  angular  deflections.  As  a 
result  of  these  successive  collisions  and  cumulative 
redirections,  one  may  want  to  compute  the  probability  that 
the  particle  emerges  from  the  slab  with  a  particular 
lateral  displacement  and  direction. 
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As  a  measure  of  multiple  scattering,  the  ICRU(17) 
suggests  the  mean  square  scattering  angle,  ©*  and  defines 
the  quotient  0VP'5L  as  the  "mass  angular  scattering 
power"  in  a  manner  analogous  to  mass  stopping  power.  This 
quantity  should  be  used  cautiously,  however.  The  original 
calculation  of  Rossi(21)  used  by  ICRU,  was  based  on  the 
validity  of  superposing  numerous  small  and  independent 
angular  deflections.  This  is  generally  not  the  case  in 
multiple  successive  scattering  processes,  especially 
through  thick  slabs.  It  is  valid  only  when,  for  each 
collision,  the  direction  of  the  emergent  electrons  does  not 
differ  significantly  from  the  initial  direction,  and  the 
electrons  remain  "bunched"  and  essentially  parallel  to  the 
incident  beam  direction.  Thus,  the  appplicabili ty  of  the 
mass  angular  scattering  power  is  confined  to  such  a  case  as 
the  passage  of  electrons  through  an  ultrathin  plate  of 
scattering  medium.  This  is  often  overlooked,  and  many 
authors  continue  to  misuse  this  quanti ty ( 38 , 40 )  for  thicker 
absorbers . 


1.3  Electron  Beam  Treatment  Planning 

The  principle  of  electron  beam  treatment  planning  is 
simple.  The  intent  is  to  optimize  the  treatment  of  a 
malignant  tumor  by  delivering  the  maximum  dose  to  the 
target  volume,  while  sparing  the  adjacent  healthy  tissues. 
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Thus,  the  treatment  planning  process  can  be  regarded  as  the 
procedure  of  getting  precise  information  of  the  patient 
anatomy  and  selecting  the  optimum  field  arrangement  to 
administer  a  prescribed  dose  distribution. 


1.3.1  Description  of  The  Patient 

The  necessary  patient-specific  data  for  treatment 
planning  are  the  patient's  external  contour,  the  internal 
anatomical  structures,  and  the  location  of  the  tumor.  From 
a  physical  standpoint,,  these  essential  data  are  available 
through  X-ray  computerized  tomography (CT ) .  It  should  be 
also  noted  that  the  CT  scanner  is  used  to  obtain  tissue 
density  maps,  although  the  potential  exists  to  also  obtain 
tissue  atomic  numbers(ll). 


1.3.2  Dose  Distribution 

Once  the  anatomical  data  for  a  particular  site  is 
acquired,  the  next  step  in  treatment  planning  is  to  select 
the  suitable  type  of  electron  beam(s)  which  will  provide 
the  optimum  irradiation.  The  accurate  computation  of  the 
dose  distribution  is  essential  and  the  method  used  must  be 
reliable  for  a  wide  range  of  beam,  and  patient 


configurations . 
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For  a  tumor  located  at  depth,  treatment  beams 
unavoidably  deliver  undesired  radiation  to  the  healthy 
tissues  overlying  the  tumor  in  the  treatment  region.  In 
this  context,  the  skin  (surface)  dose,  relative  to  the 
maximum  dose  has  important  implications  for  radiation 
protection  of  the  healthy  tissues.  The  typical  skin  or 
surface  dose  for  a  broad  beam  of  electrons  is  approximately 
85%  to  95%  of  the  maximum  dose,  and  varies  according  to  the 
method  of  producing  the  macroscopic  beam.  Electron  beams 
generated  by  scanning  a  pencil  beam,  rather  than  by 
foil-scattering,  and  requring  less  secondary  collimation 
reduce  the  skin  dose  significantly. 


1.3.3  Characteristics  of  Depth  Dose  Curves 

The  relative  depth  dose  curve  of  a  broad  beam  of 
electrons ( see  Figure  1)  can  be  divided  into  several  zones 
of  interest;  the  initial  build-up  region,  the  depth  of 
maximum  dose,  the  fall-off  of  dose,  and  the  "tail"  of  the 
curve  caused  by  photon  contamination.  All  of  these  are 
influenced  by  the  parameters  of  the  beam,  namely,  the 
nominal  energy,  the  field  size,  the  energy  spectrum,  and 
the  field  uniformity  and  symmetry.  What  they  have  in 
common  in  influencing  the  shape  of  the  depth-dose  curve  is 
the  degradation  and  redirection  of  electron  energy,  that 
is,  the  broadening  of  the  energy  and  angular  distribution 
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of  the  electron  beam. 

The  initial  build-up  region  of  the  depth  dose  curve 
is  where  the  initial  interaction  of  electrons  occur  with 
the  scattering  medium.  This  region  is  subject  to  the 
fluence  of  the  incoming  electrons  from  outside  of  the 
medium  and  some  electrons  backscattered  from  the  deeper 
regions  of  the  medium.  These  electrons  may  also  contribute 
partially  to  the  skin  dose,  if  they  manage  to  reach  the 
surface  of  the  medium.  Thus  every  electron  beam, 
regardless  of  its  size,  has  a  build-up  region,  even  though 
it  is  not  as  evident  for  a  very  narrow  beam. 

The  build-up  region  has  a  significant  meaning  in  the 
treatment  of  the  tumors  near  the  skin  surface,  such  as 
lesions  of  the  ear,  lip,  and  cheek.  In  these  cases,  it  is 
advantageous  to  employ  lower-energy  beams  but  it  is 
sometimes  still  necessary  to  protect  underlying  tissues  by 
placing  a  shield  behind  the  treatment  volume.  Certain 
studies  have  indicated  increases,  by  as  much  as  80%,  due  to 
back-scattered  electrons  from  such  lead  "shields".  This 
overdosage  can  be  minimized  by  using  a  shield  constructed 
of  lower  atomic  number  materials,  or  lining  the  lead  shield 
with  a  thin  layer  of  wax ( 22 , 23 , 24 ) .  The  following 
explanation  can  be  given  for  these  protective  effects. 
Firstly,  the  effectiveness  of  lower  atomic  number  material 
might  be  caused  by  the  reduction  in  the  fluence  of 
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backscattered  electrons,  as  confirmed  by  Baily.(25)  The 
lining  effect  might  be  explained  by  two  phenomena  happening 
in  the  lining  and  in  the  lead  shield.  The  electron  beam 
first  enters  the  wax  lining,  and  exits  with  a  high  residual 
energy  because  the  lining  is  comparatively  "transparent”. 
Once  this  beam  enters  the  lead  portion,  part  of  the 
electrons  are  backscattered  with  reduced  energy  after 
undergoing  multiple  scattering  in  the  lead.  The  lining  on 
the  incoming  face  of  the  lead  shield  is  then  more  "opaque" 
to  these  lower  energy  electrons  backscattered  from  the  lead 
and  reduce  the  dose  due  to  backscattering .  The  above 
possible  explanations  are  tentative  and  more  understanding 
of  scattering  in  general  is  required  for  validation. 

Another  significance  of  the  build-up  region  can  be 
found  in  the  distortion  which  occurs  when  the  beam  is 
incident  obliquely  on  the  surface.  In  this  case,  the  dose 
at  shallow  depths  increases  and  the  dose  at  deeper 
locations  decreases .( 26 )  The  increase  of  relative  dose  at 
shallow  depth  is  obvious;  while  most  of  the  incident 
electrons  are  retained  in  the  medium,  some  of  the 
multiply-scattered  electrons  can  escape  through  the  medium 
surface;  This  represents  a  pure  gain  for  the  build-up 
region  and  a  pure  loss  for  the  deeper  region .( Figure  4) 
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Surface  of  Medium 


Fig.  4. 


The  effect  of  tilted  surface  of 
a  medium  on  the  dose  distribution 
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When  the 

incident 

beam 

is  sufficently 

large  or 

the 

electrons 

in 

the  beam 

are 

sufficiently 

degraded , 

the 

electron 

beam  delivers 

the  maximum  relative  dose  at  a 

certain 

depth 

.  This 

depth 

of  maximum 

dose ( dmax ) 

is 

believed 

to  be 

the  depth 

where 

the  transport 

and  the 

loss 

of  electrons  are  combined  to  provide  a  maximum  net  number 
of  electrons  per  unit  volume. 

The  dose  fall-off  region  of  the  depth  dose  curve  is 
believed  to  be  caused  primarily  by  the  stopping  of  the 
incident  electrons  after  progressive  energy  loss. 
Historically,  the  implementation  of  electron  beams  for 
radiation  therapy  was  accelerated  by  the  need  for  a 
fall-off  of  dose  behind  the  target  volume.  For  lower 
energy  electron  beams,  this  dropoff  is  very  abrupt  and 
useful  for  the  treatment  of  superficial  or  shallow  lesions. 

For  monoenergetic  beams  produced  by  a  scattering 
foil  and  collimation,  the  production  of  bremsstrahlung  in 
these  devices  is  unavoidable.  The  bremsstrahlung  usually 
forms  the  tail  end  of  the  depth  dose  curve.  The  relative 
magnitude  ranges  up  to  10%  of  the  maximum  dose,  depending 
on  the  accelerator  energy  of  electrons.  Few  experimental 
studies  of  this  "tail"  and  its  significance  on  treatment 
planning  have  been  reported. 
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1.3.4  Field  Shaping 

The  treatment  of  tumors  with  an  electron  beam 
sometimes  requires  considerable  shaping  of  the  irradiation 
area.  Lead  sheets  cut  in  the  shape  of  the  treatment  portal 
are  often  used  in  conjunction  with  the  standard  circular  or 
rectangular  applicators.  The  thickness  of  lead  shielding 
is  usually  chosen  to  reduce  the  dose  in  the  "shadowed" 
region  to  less  than  5%  of  the  dose  in  the  open  beam.  For 
high  energy  electrons  additional  thickness  beyond  a  few 
millimeters  does  not  decrease  the  dose  further  because  of 
the  bremsstrahlung  produced. 

It  is  generally  suggested  that  if  the  field  is  large 
and  the  blocking  is  minimal ,  then  the  normal  dose 
calibration,  may  be  used,  but  for  extensive  or  complex 
blocking,  a  special  calibration  should  be  carried  out(33). 


1.3.5  Isodose  Curves 

In  electron  beam  treatment  planning,  isodose 
distribution  data  play  a  key  role  in  that  they  provide  the 
two-dimensional  map  of  the  dose.  Information  about  the 
dose  at  off-axis  points  and  at  the  edges  of  the  collimated 
beam  is  thereby  obtained.  It  is  common  practice  to  choose 
a  beam(s)  so  that  the  target  volume  lies  entirely  within 
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the  80%  or  90%  isodose  curve.  However,  even  in  a 
homogeneous  medium,  the  isodose  curves  for  electron  beams 
produced  by  foil-scattering  are  not  "flat”  so  that  the  area 
covered  by  80%  or  90%  isodose  line  is  considerably  less 
than  that  based  on  the  geometric  field  area  defined  at  the 
surf ace .( Figure  5)  Therefore,  to  cover  the  target  area 
completely,  a  larger  field  at  the  surface  is  usually  used, 
thereby  increasing  the  irradiation  of  healthy  tissues. 
This  problem  can  become  worse  if  tissue  inhomogeneity  is 
present.  For  example,  the  "bulging"  of  isodose  can  be 
accentuated  in  a  low  density  region. 

1.3.6  Interactions  in  Heterogeneous  Tissue 

For  the  energy  range ( <50MeV)  of  electrons  used  in 
radiation  therapy,  the  most  important  microscopic 
interactions  of  electrons  with  the  constituent  atoms  of  the 
medium  are  due  to  electromagnetic  forces.  The  number  of 
interactions  among  the  incident  "bare"  electrons  is 
negligible  compared  with  that  between  incident  electrons 
and  atoms  of  the  medium.  The  two  dominant  resultant 
effects  of  the  interactions  are  deflections  and  the  loss  of 
kinetic  energy(21).  The  deflection  takes  place  when  an 
electron  passes  in  the  neighborhood  of  a  nucleus.  Since 
the  nucleus  is  very  heavy  with  respect  to  the  incident 
electron,  the  energy  of  the  nucleus  does  not  change 
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Energy:  16.0MeV 

SSD:  100  cm 

Field  Size:  10x10  cm 


Fig.  5.  Gonstriction  of  the  80%  and  90%  isodose 
line  in  comparison  with  the  geometric 
field  size  (at  surface) (33) 
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appreciably  by  recoiling  and  the  collision  is  an  elastic 
one(exclusive  of  the  radiation  emitted  during  the  process). 
The  loss  of  kinetic  energy  occurs  when  an  incoming  electron 
interacts  with  one  of  the  atomic  electrons  having 
comparable  mass. 

From  a  physical  standpoint,  the  inhomogeneities  in 
the  irradiated  volume  can  be  regarded  as  a  distribution  map 
of  materials  having  high  or  low  density  and/or  atomic 
number.  In  case  of  human  tissues,  where  there  is  not  a  big 
difference  in  mass  numbers  of  major  constituent 
atoms (H , C , 0 , N ) ,  the  density  variations  affects  multiple 
scattering  of  electrons  by  changing  the  possible  number  of 
scattering  events,  either  elastic  or  inelastic.  That  is, 
the  density  effects  of  an  inhomogeneity  can  alter  both  the 
gross  amount  of  energy  deposition  by  increasing  or 
decreasing  the  frequency  of  inelastic  scatterings,  and  the 
overall  spatial  distribution  of  energy  deposition  by 
increasing  or  decreasing  the  number  of  elastic  scatterings. 
These  two  effects  are  interrelated  since  the  energy  loss  of 
individual  electrons  by  inelastic  events  affects  the 
elastic  scattering  which  is  energy-dependent.  The  major 
influences  of  the  atomic  number  of  the  material  is  through 
dependence  in  the  microscopic  scattering  cross-section. 
Since  the  inelastic  scattering  cross-section  contains  the 
term,  Z/A,  there  is  not  a  big  difference  in  the  collisional 
energy  loss  in  each  collision  with  an  orbital  electron  for 
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two  media  of  different  atomic  numbers.  On  the  other  hand, 

the  elastic  scattering  cross-section  contains  the  term, 

2 

Z,  /A,  in  which  variation  of  Z  can  result  in  a  significant 
effect.  The  main  influence  of  variation  in  Z,  therefore, 
is  the  angular  redirection  of  electrons.  In  summary,  the 
density  effect  of  the  inhomogeneity  alters  both  the  gross 
amount  of  energy  deposited  in  the  finite  volume  and  its 
spatial  distribution,  whereas  the  atomic  number  mainly 
alters  the  latter  only. 


1.4  Present  Theoretical  Model  for  Electron  Transport 

The  interest  in  the  transport  of  high  energy 
electrons  in  a  scattering  medium  has  been  steadily  growing 
ever  since  the  concept  of  "pencil  beam  analysis"  was 
proposed  by  Lillicrap  et.  al.(41)  Some  semi-analytic  and 
empirical  approaches  to  the  transport  problem  which  had 
been  used  previously  were  not  suitable  for  narrow  beams  of 
electrons ( 47 )  because  they  relied  on  broad  beam  parameters. 
These  methods  are  reviewd  in  chapter  II. 

A  transport  equation,  currently  used  widely  in 
medical  radiation  physics,  is  that  derived  by 
Fermi. (46)  This  equation  was  intended  for  the  description 
of  very  high  energy GeV )  charged  particles  passing 
through  a  thin  slab  of  scattering  medium.  Its  derivation 
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involved  approximations,  which  are  not  necessarily  suitable 
for  a  thicker  media  and  lower  energy  beams.  The  two 
notable  approximations  are  the  assumptions  that  particles 
(a)  do  not  lose  their  energy  in  crossing  a  medium,  and  (b) 
always  travel  paraxially  along  the  incident  direction.  The 
solutions  to  the  Fermi  equation  clearly  shows  the  effect  of 
the  latter  approximation . (46 )  The  resultant  angular  and 
spatial  distribution  functions  are  Gaussian  in  shape;  the 
probability  that  a  particle  has  large  angular  orientation 
is  negligible.  In  this  work,  the  problem  of  electron 
transport  is  approached  by  studying  a  more  general 
equation,  namely  the  linearized  Boltzmann  equation. 


CHAPTER  II 

REVIEW  OF  DOSE  CALCULATION  METHODS 


Accurate  calculation  of  the  dose  distribution  for  a 
specific  clinical  situation  is  one  of  the  main  steps  in 
electron  beam  treatment  planning.  After  the  introduction 
of  electrons  as  a  mode  of  radiotherapy,  the  dose 
calculation  methods  have  been  steadily  growing  in  number  as 
well  as  in  improved  accuracy (47 ) .  Having  started  with  the 
assumption  of  a  homogeneous  patient,  today's  dose 
calculation  methods  attempt  to  handle  the  complex 
inhomogeneity  distribution  of  tissue  in  a  patient.  The 
advent  of  X-ray  computerized  tomography ( CT )  and  the 
availability  of  high  speed  computing  facilities  stimulated 
the  development  of  more  accurate  and  yet  practical  methods 
of  dose  calculation. 

In  general,  the  early  dose  calculation  methods  were 
developed  for  application  to  simple  situations.  Current 
trends  are  towards  establishing  a  flexible  method  which  is 
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useful  for  more  complex  situation  in  heterogeneous  tissue. 
In  the  following  sections,  the  major  dose  calculation 
methods  are  reviewed (47  ) . 


II. 1  Absorption  Equivalent  Thickness  Method 

As  previously  explained  in  Chapter  I ,  the  changes  in 
density  and/or  atomic  number  of  a  scattering  medium  alter 
the  depth  dose  curve  for  a  broad  beam  of  electrons. 
Laughlin  (50)  discovered  some  general  aspects  of  these 
alterations  (Figure  6); 

(a)  Within  and  beyond  tissues  of  lower  density,  such  as 
lung,  the  central  axis  dose  increases,  as  if  the  depth 
dose  curve  is  "shifted"  longitudinally  towards  a  deeper 
location . 

(b)  Beyond  high  density  regions,  such  as  bone,  the  dose 
contours  are  displaced  toward  the  skin  surface  with  a 
consequent  decreased  dose  at  a  given  depth. 

(c)  The  magnitude  of  the  actual  dose  in  and  beyond  low 
density  regions,  relative  to  that  in  a  unit  density 
medium,  depends  in  a  nonlinear  manner  on  the  location, 
extent,  and  density  of  the  inhomogeneous  region. 
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Depth 


Fig.  6.  Central-axis  depth  dose  distribution  produced 
in  (above)  high- impact  polystyrene  and  (below) 
2.0  cm  high- impact  polystyrene  followed  by 
cork,  by  electron  beams  of  indicated  energies. 

(6) 
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(d)  Within  high  atomic  number  materials,  such  as  bone,  and 
immediately  beyond  the  interface  region  there  can  be 
regions  of  elevated  dose  ("hot  spots")* 

A  correction  factor  has  been  developed  to  take  into 
consideration  the  presence  of  inhomogeneities  and  is 
referred  to  as  the  absorption  equivalence  thickness  (AET) 
method.  Although  the  correction  factor  should  depend  on 
depth,  for  practicality,  an  average  factor  is  used  based  on 
the  following  relations; 

AET  =  1.3  Pa 

where  is  the  assumed  lung  density.  For  bone, 

.  _  Pb  (electrons/g  in  bone) 

AETbone  ~  J  - - - 

1.0  g/cm  (electrons/g  in  water) 

Pb  3.0*1023/g 
1.0  g/cm3  3.36X1023/g 

=-  0.88  Pb 

ii 

where  Pb  is  the  density  of  a  bone. 

The  correction  of  an  isodose  distribution  in  a 
homogeneous  medium  using  the  AET  factor  is  usually  made  by 
assigning  AET  factors  to  each  inhomogeneous  region  along 
divergent  rays  from  the  virtual  source.  The  corrected 
isodose  distribution  can  be  obtained  by  joining  equal  dose 
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ANTERIOR 


Fig„  7. 


Calculation  of  isodose  distribution  with, 
correction  for  bone  and  lung  inhomogeneity . 
(50) 
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II. 2  Coefficient  of  Equivalent  Thickness  (CET)  Method 


Another  way  to  determine  the  dose  distribution  in 
lung  uses  an  appropriate  correction  factor,  called 
Coefficient  of  Equivalent  Thickness (CET) ,  which  can  be 
defined  by  referring  to  Figure  8.  Considering,  first,  the 
case  where  there  is  no  divergence  of  the  beam,  the 
coefficient  of  equivalent  thickness  is  defined  as  the  ratio 
of  the  di  stance  in  the  water  to  the  distance  in  the 
inhomogeneity,  as  measured  from  the  interface,  at  which  we 
have  the  same  percent  depth  dose  value. 

CET  =  x./x* 


where 

8x(>o  =  y,  cx.) 


Here,  Q  refers  to  percent  depth  dose,  X  to 
measured  from  the  interface,  and  1  and 
standard  water  curve  and  water-cork  curve, 
These  have  to  be  measured  experimentally. 


the  distance 
2  refer  to  the 
respect i vel y . 


Now  the  decreasing  part  of  the  curve  in  water  can  be 
described  by  the  empirical  equa t ion ( 5  1 )  ; 

y  =  1 1  o  -  io  eHX 

where  x  is  measured  from  the  depth  of  the  100%  point  on  the 
central-axis  depth  dose  curve.  The  value  can  be 
calculated  from  the  standard  water  curve; 


' 
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DEPTH  (cm) 


X 


Fig.  8-  Coordinates  and  parameters  used  for  the 
calculation  of  GET  and  MAC  methods.  (47) 
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y,(x,)  -  no- 10 e 


H(*i  +  O 


where  c  is  the  distance  between  the  interface  and  the  depth 
of  the  100%  point(dmax)  on  the  water  curve;  c  can  be 
either  positive  or  negative; 


(Xa-CET  +-  C  ) 


For  a  given  ray,  there  will  be  an  inverse  square  correction 
given  by 


5S&  +  d.  1-CE-T-Xx 

SSD+  do  +•  Xi 


and  Hitxp  =  y  ,(Xi)  '  S 


Thus,  the  overall  expression  is; 

|4  (X./CET  +  C)l(  SSP  +  -t-CET-Xx 

f  \  SSD  +cie  +  Xi 


Since  is  measured  at  Xz  >  H  can  be  found  from 
the  water  depth  dose  curve,  and  since  SSD,  ,and  c  are 
known,  the  equation  can  be  solved  for  CET. 

As  in  the  AET  method  described  previously,  the  beam 
can  be  divided  into  seperate  rays,  a  new  attenuation 
coefficient  |A  and  values  for  x  and  d6  measured,  and  the 
corrected  doses  calculated  along  each  ray. 
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II. 3  Modified  Absorption  Coefficient  (MAC)  Method 

To  implement  the  CET  method  with  consideration  for 
the  polarization  effect,  field  size  dependence,  deviation 
of  central  axis  percent  depth  data  from  Laughlin’s 
equation(51)  at  electron  energies  above  20  MeV ,  and  the 
variation  in  the  densities  of  internal  organs  for  the 
specific  patient,  Bagne(52)  introduced  the  Modified 
Absorption  Coef f icient(MAC)  Method. 

For  a  given  percentage  depth  dose  measured  in  a 
water-cork  phantom  at  a  distance  Xz  from  the  interface,  the 
MAC  value  is  defined  as  (see  Figure  8); 

A\AC  =  x^/XzA 

where 

and 

5SP  +  d,  +MAC-Xzpz/p, 

SS  D  +  d0  +  Xj. 

Here,  the  water  and  cork  densities  are  designated 
by  Pj  and  ,  repectively,  and  S  is  the  inverse  square 


correction  factor. 


t 
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Although  percentage  depth  dose  curves  for 
energies  <  20  MeV  follow  Laughlin's  equation,  data  obtained 
for  higher  energies  deviate  significantly  from  a  straight 
line  at  values  below  the  50%  percent  depth  dose  level(52). 
Bagne  suggested  a  modified  equation  based  on  curve  fitting; 

I (x) =  110-10  e*p{p[(x-dm)P  -  \(X-dmy  P5]  }  S 

where  p  is  the  density  of  the  medium  and  X  is  the  depth 
dose  correction  coefficient.  Regardless  of  field 
size,  X  is  constant  for  a  given  material  and  energy.  For 
energy  20  MeV,  \  approaches  zero. 


The  percentage  depth  dose  at  a  given  point  in  water 
is  then  obtained  as; 


=  (it o  -  io<2^p{o. P  , 


where  P  is  the  correction  for  the  polarization  effect,  c  is 
the  distance  between  the  100%  depth  dose  point  and  the 
surface-interface  distance,  and  n  is  defined  as; 

fl  ~  (  X  ~  d  m)  /  (  Sc\0  ~  dm) 


Then,  yz(Yz)  can  be  expressed  in  the  following  way  using 
the  MAC  method; 

l(Xa)  =  {  110  -  10e*p[ 0.M3  ED/(d,0-cU)]}  SP , 


ED  =  (MAC  •  CP,)  -X(MAC-P.-Xr-  CP,)? 


where 
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I *1.4  Correction  Method  for  Edge  Effects 


For  the  determination  of  the  dose  perturbation  at 
the  edge  of  a  finite  inhomogeneity,  Pohlit(54)  introduced  a 
correction  method  based  on  experimental  observations »  In 
order  to  quantify  the  limit  of  influence  of  the  edge, 
Pohlit  graphically  defined  two  angles  ol  and  Q  ,  which  are 
dependent  on  the  mean  electron  energy  E  (Figure  9  and 
Figure  11 ) ; 

E  =  Eo  (I-  2/Rp)  , 

where 

—  initial  energy , 

=  practical  range, 

2  =  depth  in  medium. 


The  maximum  change,  Pm^X  ,  of  absorbed  dose  due  to 
an  edge  between  different  kinds  of  medium  is  (Figure  10); 


where  t'j'ttixx  is  the  absorbed  dose  at  the  highest  increase  or 
depression  due  to  the  presence  of  the  edge  and  P0  is  the 


absorbed  dose 

at  the 

same  point 

in 

the 

absence 

of 

the 

edge . ( Figure 

11) 

Whenever  edges 

are 

present 

in 

the 

irradiated  volume ,  the  dose  can 

be 

corrected  along 

the 

angles  OL  and  0  . 


. 
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Fig.  9.  Energy  dependence  of  the  angles  ct and  £. 
(54) 


Fig.  10.  Maximum  changes  in  absorbed  dose 
behind  edges.  (54) 
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Fig.  11.  Measured  dose  distribution  behind  a  lead 
strip  in  a  water.  (54) 
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II. 5  Approximation  of  Isodose  Curves  by  the  Age  Diffusion 
Solution 

For  the  accurate  characterization  of  the  dose 
distribution  of  rectangular  electron  beams  with  minimal 
amount  of  empirical  data,  the  solution  to  the  age-diffusion 
equation  can  be  used  to  approximate  the  dose  distribution, 
at  least  in  a  homogeneous  medium(36). 


One  typical  example  is(see  Figure  12); 


}] 


*  f/(f+h)J  x  ^f[-(2y?Rp)-(fcT)'A 


where 


D  (Xj3,2,T)  =  dose  at  depth  Z, 


p  = source- skin  distance. 


X0  —  (  1/2  )  (  Wi  cH"h  4- Ke  )  •  £(  F  h- 2)/Fj  ? 

go  -  0/z)  (  l-eng-FF+  Ke)  °[(F+2)/f]  , 


INp  J  ' 

Q,  ,  C,  M,  p  =  constant  for  a  given  energy. 


D  —  practical  range. 


' 


' 
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Central  Axis(cm) 


(a) 


Central  Axis(cm) 


6  A  2-  0  2  4  $ 


(b) 


Fig.  12.  Calculated  isodose  distributions  in  water, 
assuming  a  monoenergetic  incident  electron 
energy  of  15  MeV  and  a  field  size  of  10  cm 
with  a  parallel  beam.  The  values  of  c  and 
n  are:  (a)  Left  side,  c=1.5,n=1.0  Right 
side,  c=1.5,n=3.0  (b)  Left  side,  c=1.2,n= 
2.0.  Right  side,  c=1.8,n=2.0.  (36) 
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All  the  contants  are  to  be  "adjusted"  to  give  a  good 
fit  to  the  experimental  data.  The  functions  of  individual 
parametric  constants  are;  P  characterizes  the  beam 
diffuseness  at  the  surface;  C  characterizes  the  increase 
in  diffuseness  with  depth;  N  is  an  overall  shape-changing 
parameter;  Ke  adjusts  the  width  and  length  of  rectangular 
beams.  This  method  is  used  on  some  commercial  treatment 
planning  computers (e . g . ,  AECL  TP-11),  but  no  allowence  are 
made  for  tissue  inhomogeneity. 


II. 6  Pencil  Beam  Analysis  of  Broad  Beam 


Unlike  most  of  the  dose  calculation  methods, 
discussed  previously,  pencil  beam  analysis  method  deals 
with  a  broad  beam  of  electrons  not  as  a  whole  but  as  an 
integration  of  a  number  of  narrow  beams  (Figure  13).  In 
this  way,  dose  calculation  problem  becomes  a  study  of  a 
narrow  beam  of  electron  which  is  simpler  in  many  aspects 
(e.g.,  mathematical  treatment)  and  can  be  used  to  account 
for  small  inhomogeneities. 


As  an  analytic  expression  for  a  narrow  beam,  authors 
have  used  a  solution  to  the  Fermi  equation ( 55 , 56 ) „ 

eHz)  x  -  2.  re  (z)x  e*  +-  r(z)  e*  1 


{(Qx,X,i)  = 


©X.F> 


©l(z')  roo  - 1  ve(z) 


TT  {_©*(*)  roo-  [reno]1} 


i/z 
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Z  -0 


X 


w 


Fig.  13.  Schematic  representation  in  the  X-Z  plane 
of  a  therapeutic  electron  beam  incident 
on  a  patient.  (38) 
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where 

_____  _____  __  /  i 

rzd)  -  r(o)+2reti)z  +  e"(o)zi+  \  (2-ufT(u)du 


and 

re(i)  =  re  io)  +  e?(o}z  +  (  (z-u)T(codu 

'  o 


With  this  solution  as  a  basis,  a  broad  beam  of 
electrons  can  be  represented  mathematically  as  a 
convolution ; 

<£(y,y,*)  =  f(  A (*', y', z')-f (x-x, y-y, s'-z) dx'dy ' 

J^s 

where  Au'.y1  is  the  incident  beam  fluence  and  S  is 

the  cross-sectional  area  of  the  broad  beam. 


I I . 7  Summary 

As  mentioned  earlier,  the  earlier  dose  calculation 
methods,  such  as,  AET,  CET,  MAC  mehtods,  were  developed  for 
a  specific  treatment  situation  in  which  laterally  broad 
inhomogeneities  are  present  in  the  treatment  volume(e.g., 
lung).  Accordingly,  the  applicability  of  those  methods  to 
a  more  complex  si tuation (e . g . ,  small  voxels  of 
inhomogeneity)  remains  questionable,  even  though  some 
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institutions  accept  these  methods  to  develop  a 
"comprehensive"  electron  beam  treatment  planning 
system( 52 , 55 ) .  Furthermore,  the  empirical  nature  of  the 
approaches  to  the  problem  of  these  methods  makes  their 
application  to  a  "beam  model"  for  computerized  treatment 
planning  difficult.  This  trend  is  also  true  for  Pohlit's 
correction  method  for  the  edge  of  an  inhomogeneity  and 
Kawachi's  isodose  calculation  method. 

Considering  these  limitations,  the  pencil  beam 
analysis  of  broad  beam  has  two  distinctive  advantages  for 
the  computerization  of  electron  beam  treatment  planning. 
First,  it  would  require  only  a  minimum  amount  of  empirical 
data  for  characterization  of  the  electron  beam  from  a 
particular  medical  accelerator.  Secondly,  as  the  broad 
beam  is  considered  to  be  made  up  of  a  number  of  narrew 
beams,  the  treatment  of  inhomogeneities  can  be  done  on  a 
voxel  by  voxel  basis,  using  CT  data.  However,  at  current 
stage,  the  necessary  theoretical  model  for  the  transport  of 
a  narrow  beam,  even  in  a  homogeneous  medium,  is  not  well 
established.  The  only  present  model  we  have  is  Fermi's 
model  which  is  quoted  in  section  II. 7.  This  model,  for 
example,  ignores  the  existence  of  backscattered 
electrons ( see  Chapter  II).  Furthermore,  before  we  can  use 
the  pencil  beam  analysis  as  a  "dose"  calculation  method, 
the  relationship  between  the  electron  fluence  distribution 
and  the  resultant  dose  distribution  should  be  well 
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established,  because  theory  only  provides  the 
fluence  distribution. 


electron 


CHAPTER  III 

TRANSPORT  OF  ELECTRONS  IN  A  SCATTERING  MEDIUM 

III.l  Introduction 

The  physical  problems  in  electron  radiation  therapy 
are  largely  due  to  the  beam  production  techniques  commonly 
used.  For  instance,  the  scattering  foil,  beam-defining 
collimator,  and  added  shielding  materials  are  used  to  shape 
a  field  but  these  materials  degrade  the  beam.  This  is  at 
the  origin  of  two  of  the  most  important  current  problems, 
that  is,  tissue  inhomogeneity  and  dose  standardization. 

From  a  physics  standpoint,  the  tissue 
inhomogeneities  are  local  clusters  of  material  which  have  a 
mass  density  and/or  atomic  number  different  from  that  of 
water.  These  clusters  are  usually  small  in  size  compared 
to  the  electron  beam  dimension,  and  yet  are  not  negligible 
in  size.  The  absorbed  dose  depends  on  the  relative 
location  of  the  inhomogeneities,  with  respect  to  a  point  of 
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interest.  The  assessment  of  the  influence  of  the 
inhomogeneities  on  a  broad  beam  of  electrons  is  very 
difficult.  Even  after  assuming  that  the  required 
perturbation  effect  of  the  inhomogeneity  is  characterized 
empirically,  current  schemes  of  electron  dose  computation 
do  not  have  the  necessary  foundation  to  predict  these 
effects  and  provide  guidance  for  dose  compensation. 
Currently  compensation  for  inhomogeneities  is  only  achieved 
approximately  by  selecting  an  appropriate  electron  beam 
energy  at  the  cost  of  other  advantage (e . g . ,  skin-sparing). 

The  problem  of  standardization  is  as  serious  as  the 
problem  of  inhomogeneity.  Without  standardization  there  is 
restricted  sharing  of  physical  or  clinical  results  because 
measurements  taken  on  one  medical  accelerator  with  the 
"same  energy"  may  not  be  relevant  to  another.  Beam 
spreading  by  a  scattering  foil  was  developed  without 


serious  consideration 

of  electron  transport 

and 

makes 

it 

very  difficult  to 

develope  shareable 

data 

bases 

of 

"standard  electron  beams". 

Recently,  several  attempts  were  made  to  understand 
the  characteristics  of  the  dose  distribution  of  a  broad 
beam  of  electrons  by  decomposing  it  into  a  number  of  narrow 
beams . ( 40 , 41 )  This  type  of  analysis  is  advantageous  as  a 
method  of  understanding  the  inhomogeneity  effect.  The 
lateral  size  of  a  pencil  beam  is  virtually  zero  at  the 
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surface  of  a  medium,  and  does  not  increase  significantly  at 
depth  in  comparison  with  the  size  of  inhomogeneities. 
Therefore,  assuming  the  synthesis  of  a  broad  beam  from  a 
finite  number  of  pencil  beams  is  valid,  we  need  not  be 
concerned  with  the  lateral  dimension  of  the  inhomogeneity; 
every  inhomogeneity  can  be  seen  to  be  "slab-like",  relative 
to  an  incident  pencil  beam  size. 

Pencil  beam  analysis  can  also  be  applied  to  minimize 
the  problem  of  standardization.  The  common  origin  of  this 
problem  is  the  mechanism  of  producing  a  broad  beam  of 
electrons  which  degrades  the  electrons  to  such  an  extent 
that  the  characterization  of  the  broad  beam  becomes 
ambiguous.  However,  the  characterization  of 
monodirectional  and  monoenerget ic  pencil  beams  which 
compose  the  degraded  broad  beam  should  be  possible  in  a 
unique  way,  because  the  degree  of  degradation  of  electrons 
within  each  pencil  beam  would  be  independent  and  not  too 
severe . 


III. 2  Some  Aspects  of  Electron  Transport  in  a  Scattering 
Medium 

All  the  electrons  in  a  therapeutic  electron  beam 
lose  their  initial  energy  in  a  stopping  medium,  namely,  the 
tissue  of  a  patient.  While  losing  energy  through  numerous 
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inelastic  collisions,  the  electrons  are  also  deflected  by 
elastic  collisions  with  the  nucleus.  These  two  events  are 
interdependent  because  of  progressive  degradation  in 
electron  energy  and  angular  distribution  of  the  electron 
fluence  in  a  scattering  medium.  In  the  following  section 
we  will  discuss  some  features  of  these  combined  effects. 


III. 3  Multiple  Scattering  and  the  Linearized  Boltzmann 
Equation 

Basically,  the  electron  transport  problem  in  medical 
radiation  physics  is  that  of  an  infinite  absorber.  It  is 
assumed  that  electrons,  in  general,  participate  in  a  large 
number  of  collisions  before  they  are  finally  stopped.  In 
this  calculation  we  assume; 

(a)  statistical  fluctuations (" straggling "  )  in  electron 
range  and  energy-loss  are  ignored. 

(b)  all  the  collisions  between  electrons  and  atoms  are 
assumed  to  occur  instantaneously  and  elastically. 
However,  the  mean  energy  loss  over  a  finite  path-length 
will  be  considered .( see  section  II. 3. 3.) 

(c)  bremsstrahlung  production  by  elastic  collisions  is 
ignored . 


If  we  rely  on  the  conservation  of  electrons  within  an 
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inf initisimal  volume  element  of  phase  space,  then  the 
linearized  form  of  the  Boltzmann  equation  is; 


2i 


(x.l) 


of 


where  f  &)c)a  <J\T  is  the  expected  number 

electrons  in  the  volume  element  d?  about  the  point 

j  /V  /v 

whose  velocities  lie  in  the  range  of  diT  about  AX  ,  at 
pathlength  JL  .  Y\  is  the  number  of  scattering  centers  per 
unit  volume  of  medium.  (Rtj-v,  l)  is  the  differential 


scattering  cross-section.  Here  the  scattering  centers  are 
assumed  fixed  in  position,  spherical,  and  uniformly- 
distributed  in  the  medium.  Thus,  elastic  scattering 


cross-section  only  depends  on  the  angle  between  the 
trajectories  of  incident  and  emergent 

—  |  t  A  A  j\ 

electrons ,  OSU)  VV'VJ  .  The  general  solution  to 
equation ( 1 1 . 1 )  is  uniquely  determined  if  proper  boundary 
conditions  are  specified.  The  existence  theorem  of  a 
solution  in  a  very  general  case  and  many  more  useful 
properties  of  equation  (II.  1),  have  been  studied  by  Case 
and  Zweifel.(42)  Since  the  equation  (II. 1)  is  linear,  it 
is  clear  that  the  solution  for  a  broad  surface  boundary 
condition  can  be  established  by  the  superposition  of 
solutions  with  "point"  boundary  conditions  on  that  surface. 
In  our  study,  we  are  interested  in  the  solution  to 
equat ion ( I I . 1 )  with  the  pencil  or  point  boundary  condition 
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on  the  surface  of  serai-infinite  homogeneous  medium.  In  our 
discussion,  we  will  assume  that  the  pencil  beam  is 
monoenerget ic ,  monodirectional ,  and  incident 

perpendicularly  on  the  medium  surface.  In  this  geometry, 
the  cylindrical  coordinate  system  (  P  ,  "2  )  without  an 
azimuthal  angle  for  spatial  coordinates  is  suitable.  For 
velocity  direction  coordinates,  AT  ,  spherical  coordinates 
will  be  used  throughout  discussion.  The  azimuthal  angle  of 

velocity  coordinates  is  defined  as  the  angle  between 

n  /K 

the  p  ~  l  plane  and  the  projection  of  4J  on  the  Z-plane. 

The  incident  direction  is  chosen  along  the  Z-axis. 


III. 3.1  Fokker-Plank  Approximation  At  Shallow  Depth 

—i  /  S*  /N  >  \ 

Supposing  that  the  scattering  angle  ITT*  l\Jy  is 

very  small,  the  collision  integral  within  equation ( I I . 1 ) 
can  be  approximated  by  expanding  -f  (?,  v.  50  in  a  Taylor 

A  f  ^ 

series  about  V  ,  retaining  only  the  lowest 

nonvanishing  terms.  If  we  substitute  Cf)0  0  =  |A  ,  we  find 


where  the  "diffusion  length",  D,  is  defined  by 


2TTK 
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1 

D 
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where  (tHl). 


In  order  to  estimate  the  error  involved  in  this 
approximation,  we  consider  an  electron  perpendicularly 
incident  on  the  surface  of  a  medium.  For  simplicity  we 
concentrate  on  the  probability  that  the  velocity  of  this 
electron  has  a  polar  angle  0  at  pathlength  j<L  without 
regard  to  the  spatial  coordinates.  Throughout  the  rest  of 
discussion  we  shall  be  concerned  only  with  path 
length  il  small  compared  to  the  diffusion  length  D  of  the 
medium.  For  8.«  D  ,  it  is  expected  that  f  (6,51)  will  be 
peaked  at  small  angle  0  ;  and  in  the  small  scattering 
angle  approximation,  the  equation  (II. 2)  is  solved  by 

fF-P  =  TW  ^PC-01/^)  ^ 

where 

&  =  41/D  «  i  .  UM) 

It  is  assumed  that  the  energy  of  the  particle  does  not 
change  too  much  during  the  collisions  if  Jk<D  .For  a 
simple  assessment  of  error,  we  compare  solution  (II. 3)  to 
the  true  solution  of  the  Boltzmann  equation,  using  an 
iterative  test.  By  inserting  (II. 3)  into  the  collision 
integral  of  the  Boltzmann  equation,  we  find,  with  the  small 
angle  approximation, 
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£/*p  -(^-ze^c^cpy  0^ 


(3L.5) 


In  the  Fokker-Plank  approximation  to  the  collision 
integral,  the  above  equation  reduces  to; 

TT 

dt  <rw  (i.  b) 

We  can  say  that  the  Fokker-Plank  approximation  is 
appropriate,  for  a  given  &  and  ,  if  (II. 6)  adequately 
approximates  (II. 5). 

The  equation  (II. 2)  can  obviously  be  meaningful  only 
if  the  scattering  cross-section  falls  off  sufficiently 
abruptly  with  increasing  ^  .  In  the  present  discussion, 

we  restrict  our  considerations  to  scattering  cross-sections 
which  fall  off  monotonically  with  increasing  H'*  and  which 
have  the  property  that  there  exists  some  critical  angle 
such  that 

\  (T(MO  «;  1  ^  Y  (TWO  (Tt.n  ) 

Under  this  condition,  it  is  clear  that  one  necessary 
condition  for  validity  of  the  Fokker-Plank  approximation  is 


, 
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0*  =  4S./p  >  (m.S) 

That  is,  the  approximation  is  dependable  only  after  a 
sufficient  time  has  elapsed.  Physically  this  means  each 
electron  must  have  suffered  enough  collisions  before  they 
can  be  effectively  described  by  the  equation  (II. 2).  After 
the  first  mean  collision  time  interval,  the  true 
distribution  should  already  be  "filled  out"  to  a  value  of 
order  •  Thus,  enough  time  is  needed  for  C  /  to 
exceed  »  that  is,  for  the  orientations  to  be  relatively 
random  or  "diffuse".  It  should  be  also  noticed  that  the 
expansion  of  equation  (II. 5)  to  terms  linear  in  'V  is 
incorrect  unless 

0  «  01/  fr-'O 

which  is  a  necessary  condition  for  the  Fokker-Plank 
approximation  to  be  adequate.  This  condition  too  has  a 
physical  significance.  A  typical  collision  involves  a 
scattering  through  angle  of  order  *  In  a  sequence 
of  M  collisions,  with  random  changes  in  azimuthal  angle, 
the  mean  net  deflection  angle  ( 0Z)  is  given 
by  .  On  the  other  hand,  the  improbable  sequence 
of  such  collisions  in  which  the  velocity  vector  remains 
essentially  in  a  fixed  plane  yields  a  net  angular 
deflection  e  ^  nj^o  .  Such  a  coherent  chain  of 


scatterings  is  not  properly  accounted  for  in  Fokker-Plank 
approximation.  We  therefore  expect  the  solution  to 


I 
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equation  (II. 2)  to  become  invalid 
when  0  y  M  ^  &~  /  .  While  conditions  (11.8,9)  are 
necessary  for  securing  the  validity  of  the  Fokker-Plank 
approximation,  the  sufficiency  is  provided  by  a  simple 
obvious  condition; 


e"<  e1 


(X.  ID) 


as  one  expects  on  physical  grounds  and  can  verify  by 
inspection  of  equation  (II. 6).  For  an  improved  estimate  of 
this  sufficient  condition,  we  inquire  about  the 
contribution  to  -p(0yjl)  for  ©2^  0*  .  This  contribution 
comes  from  a  single  large  angle  scattering  having  taken 
place  during  the  path  length  interval  X  .  Since  the  major 
contribution  to  the  probability  distribution  comes  from  the 
deflection  angles  up  to  a  value  of  order  0^0  ,  such  a 
scattering  must  involve  an  angular  deflection  which  is 
essentially  equal  to  0  itself.  Hence,  this  single 
scattering  contribution  is  approximately  given  by 

■fjle.ii)  cz  2.xn£cr(e)  ,  e»-e  (x.n) 


Note,  from  (II. 4)  and  (II. 11),  that 

-fs  /  -  (2Ttn/4)  DCT(e)  (X.iz) 


f 


and  from  the  definition  of  P 


that 


Now ,  from  the  single  large 
can  calculate  the  total 
scattering  through  angles 
integration  over  angle  © 
supposing  that  eVo) 
function,  we  have 


angle  scattering  probability,  we 
probability  Pt  (©.>&)  for  single 
of  order  0  or  larger  by  simple 
.  Since  ^V©^©^  0^"^  ,  and 

is  a  monotonically  decreasing 


PtO.JO  = 


/-w 


n  e‘  d  e1cne')  £ 


eVie) 

ex  iTm.) 


(n .  (4) 


Therefore,  we  can  safely  neglect  multiple  large  angle 
scattering  for  all  6»e  .  However,  single  large  angle 
scattering  can  be  neglected  only 

when  fs  l&A)  «  fp-p  .  The  equation  (II. 2)  is  then 

adequate  provided  ©  satisfies  (II. 9)  and  is  also  less  than 
a  critical  angle,  0C  ,  determined 

by  ff-r  (aJL)<-Ts  (fyl)  ■  Thus, 

—hr  &<P  (-Q^/e2)  ~  2irn  8.  (T (,0c)  (X.  15) 

IT 


or  from  (II.4)and  (11.13); 

e*p  (- of/ e1)  ~  TT  (  HC/  0c) 


(x, it ) 


. 
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These  results  can  also  be  understood  from  a  direct 
comparison  of  equations  (II. 5)  and  (II. 6).  For  0  9  , 
the  expression  in  curly  brackets  in  (II. 5),  regarded  as  a 
function  of  ,  has  a  sharp  peak  of 
height  ««*p(ey  e*)  and  of  width  0  ,  centered 
at  In  the  integral  on  the  right  hand  side  of 
equation  (II. 6)  this  makes  a  contribution  which  is  of  the 
order  of  ;  i.e.,  it  is  just  the  contribution  of 
a  single  large  angle  scattering.  Indeed,  if  the  condition 
(II. 9)  is  satisfied,  the  right  side  of  equation  (II. 5)  is 
essentially  the  sum  of  two  terms; 


2£ 

Til 


HI  H4 


In  summary;  (1)  The  transition  from  the  Boltzmann 
equation  to  equation  (II. 2)  is  acceptable  only  if  the 
scattering  cross-section  (TWO  falls  off  sufficiently 
rapidly  with  increasing  Vp  (faster  than  'p  4  )  beyond  some 
small  angle  .  (2)  Even  then,  any  solution  of  equation 
(II. 2)  must  be  thought  of  as  representing  a  smoothing  of 
the  true  distribution  function  over  path  length  intervals 
large  compared  to  Where  the  distribution 
function  has  already  become  smoothly-varying  in  phase 
space,  this  limitation  is  no  longer  important.  However,  as 
in  the  example  we  have  analyzed  the  initial  *f“  is  sharply 
peaked  in  phase  space,  the  solution  cannot  be  trusted  until 
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an  interval  of  time 

zJv 

has 

elapsed.  In 

any 

case  , 

if  (TOP)  is  slowly 

varying 

in 

the  interval 

changes  in  the 

distribution 

function 

over 

an 

interval  A/l'xJlo  can 

always 

be 

obtained  from 

pure  single 

scattering  analysis. 

(3)  If 

at 

some  initial 

time 

the 

distribution  function  is  essentially  confined  to  a  narrow 
region  in  phase  space,  its  subsequent  development  within  a 
larger  domain  of  phase  space,  over  an 
interval  X.C<i«D  ,  is  adequately  approximated  by 
regarding  as  a  direct  sum  of  the  main  solution  and  a 
contribution  “Fs  coming  from  single,  large  angle 
scatterings ; 

f  ~  fp-p  +  is  Ot.'S) 

where  "Ps  is  to  be  included  only  for  angles  large  compared 
to  the  mean  angle  for  ff~P  In  the  same  problem 
analyzed  above,  the  domain  of  validity  of  (11.18) 
corresponds  to  0  <C  0  <^C  B2-/  ,  the  term  being  appended 
only  for  0»e  and  becoming  important  only  beyond  some 
critical  angle,  @c  determined  by  setting  pp.p  ^  P5  and 
taking  the  root,  if  any,  which  lies  well  above  0  . 
Finally  a  proper  assessment  of  the  errors  involved  in  any 
approximation  scheme  for  treating  multiple  scattering  must 
depend  on  the  details  of  the  problem  in  question. 


. 


62 


III. 3. 2  Applications  to  the  Shallow  Depth  Problem 


At  the  initial  time,  an  electron  is  assumed  to  be 
located  at  the  spatial  origin  7=o  with  velocity  vector 
directed  along  the  positive  Z  axis.  In  the  following 
treatment,  we  scale  the  distances  by  measuring  them  in 
units  of  ,  the  diffusion  length.  D  is  assumed  to  be 
independent  of  in  the  region  under  consideration.  We 
start  with  studying  the  dependence  on  orientation  of  the 
velocity  vector  of  the  low  spatial  moments  of  f  ,  in 
particular  for  times  which  are  small  compared  to  the 
diffusion  time,  i.e.,  for  x«i . 


Consider  first  the  zero-th  spatial  moment; 


-f„  (0-0 


-f(P,2,e, <p.*)  NPdz 


(in) 


which  describes  the  distribution  in  velocity 
orientation,  irrespective  of  spatial  position, 
moment,  after  averaging  over  ,  equation  (II. 2) 
to 


vector 
For  this 
reduces 


7>H  r  1%  ' 


(x,20) 


and  the  solution  consistent  with  the  boundary  condition  is 

l  °9  ,  . . v  n  /|A 


(31.21) 
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where  Pm(F)  is  the  Legendre  polynomial . 


In  order  to  find  out  the  rapidly  convergent  solution 
for  51  ^  I  ,  we  exploit  the  integral  representation  of 
Legendre  polynomial;  i.e.,  for  ^  } 


^  ^  ^  )9  (cwe-cw«)'/J 


so  that 


4'TTfe  = 


TT 


dc< 


IM 


Jzir  )  ”  (aw  0  -  cXJ  (X  )!/z 


where 


CoO  =  ^T  Cim+D-<wn[^+s)c<]  e 


-m  (tyH-01 


m-  -cm 


With  the  Poisson  sum  formula,  ^.(^) can  be  transformed 


into 

S.(pO  = 


£Xi 


(X> 


,  -iTTmx  ,  r,  .  - 

axe  (2X+1)  Ain  I  (x  +  i)  <*  j  € 


X(X+02 


ex> 


~  Ml  ,  -2TrimX  -Qx2  1/4 

I  <1*  g  CM«e  e 


^  ih  y 

loc  rn  ;_,oo 


For  sufficiently  small  A  ,  only  the  1 OH— 0  term 


retained  in  the  above  expression  and  we  find 

l  nr.  c(e-0tV+JL 


4t£ 


- 


<Ef  A* 


e  (um  e  -  ok>  o(),/z 


need  be 
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We  evalute  this  in  the  limit  ^  O  (but  for  arbitrary  &  ) 
and  we  get 

I 


fo  (Q'A) 


4^5.  v  Aw\0 


9  V/z  -eV4i  ,  . 

£  (x  2?) 


which  is  consistent  with  the  earlier  result  (II. 3). 

Now,  to  get  information  about  average  longitudinal 
displacement,  we  consider  the  higher  moments  of  the 
distribution  in  the  coordinate  of  spatial  position, 
regarded  as  functions  of  the  velocity  vector  orientation. 

f  lnJn  (e,iO  s  j  z"f  (P,  2,  e,  9,  fi)  (x.  24) 


Multiplying  the  equation  (II. 2)  by  2  , integrating 
over  all  space,  averaging  over  azimuth  angle  of  velocity 
vector  orientation,  we  find  a  set  of  coupled  equations; 


Df(nji) 

Ti 


DjA 


(\-U*)  2£cn;i))  =  uH-fcn-i;?) 

I  ^  J 


(3L.25) 


This  is  to  be  solved  subject  to  the  initial 
condition  ^  0  at  a-o ,  provided  H^>0.  It  is 
evident  that  the  solution  for  ,  in  terms 
of  “fcn-l;r>  »  is  given  by 


(e,  *)  =  ;&(  d2'^(2m+-i)Pm(p)  e 
0 


•mtw+n  tx-Ji1) 
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x  Jil'  Pvti ( \x")  (b',S.[) 


(iL,  lb) 


Since  the  zeroth  moment  -p0  is  already  known,  all  of  the 
higher  moments  can  be  obtained  in  succession.  For  example, 


r  I  v  d  Mf  -mcm-OA  -Crn-nXm-t-DJO  ,  . 

=  ^  z  Pm(p)  [e  -e  ^  (1L'2n) 


i 


Using  the  Poisson  sum  formula  as  before  we  find  for  the 
limit  I  0  the  result; 


fo; 


i) 


->  A 


~e~ 


f0u,n 


Now 


Z  -  f  (Hz')  (®>  #  *0 


(1,2?) 


is  the  average  value  of  the  2  displacement  at 
"time"  iL  for  electrons  whose  velocity  vector  makes  an 
angle  0  with  respect  to  the  Z  axis;  so  we  have 


2 


©c 


-mcm-oA  -cmH)crvH-2)£ 


g,PmW{e.  -e 


CO 


(2UH-I)  Pm(fA)  g 


-mCm-H)  A. 


(1,2^) 


and 


£  0 
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We  consider  next  the  question  of  the  lateral  spread 
with  time  of  a  group  of  electrons  all  located  initially  at 
the  origin,  with  the  velocity  vectors  directed  along 
the  2  axis.  Define  the  moments 


The 


cn;  p) 


P  -f  l)  PdPd2 


/ 

equation  (II.  2)  then  reduces  to 

if  imp)  f  JL  (i-u1-)  ^fcmp)  +  _J _ 

D8.  (Tp  r  I  -jjO  . 

—  Avyi  0  C&A  ip-fui-l;  f>) 


(X.30) 


(11.30 


which,  for  r\>o 


must 


be 


solved 


subject 


to  -ftp*  p)  ~  0  as  O  .  This  has  the  solution 


•fen; 


f  A  co  \ 


P ) 


-  n 


da'r  -  v  I  Xm 

a~o  m=-\ 


Y>Je^  e 


xa+ou-i') 


X 


dY  Yx*  (eW'^ecWW P)  W. O 

(ir.  32) 


where  the  (Xm  are  spherical  harmonics.  To  determine 
moment  functions  we  proceed  as  before.  With 
given  “fo  (Bj  9-)  ,  we  find 

0°  *  -V 

Pw(e)] 


the 

the 


Yup)  (©.‘P’  5)=  X 


l 


£7T  10  p,=,  m 
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r  -wum+oa. 

Me  -  e 


(U.  33) 


Thus  , 

P(e,<P,Jl) 


fo  ( &,  s.) 

3  1 

OMCp-jQ  ^  rn 


{p„-1(8)-f„(e.)}{e’'“7*e 


mcm+iX" 

J 


zZ  (zm+Opmle)  e'^04- 

(tl34-") 

To  get  the  behavior  of  p  for  the  limit,  ,  we  work 

directly  with  the  differential  equation  (11.31).  We  set 


f(,;p>  (e,cp,«.)  =  P(e.Mfi)-foCe»si) 


(X  ■  35) 


and  insert  this  into  equation  (11.31)  and,  in  the  resulting 
equation  for  P  ,  retaining  only  the  linear  term  in  an 
expansion  in  power  of  5L  .  The  result  is  given  by; 


P (©,?,  2) 


0 


©►  SloskJcpi 


i-  owe 

© 


(H3b) 


the  interpretation  of  this  is  the  following.  Let  P*  be  the 

projection  on  X  -  ^  plane  of  the  position  vector  ;  and 

— 

let  irf  be  a  unit  vector  directed  along  the  projection 
ys. 

of  “IT  in  the  X  -  ^  plane.  Then  the  average  value  of  p  is 
expressed  by 
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III. 3. 3  Diffusion  In  Deeper  Regions 


In  the  deeper  region  of  a  medium,  the  variation  of 
elastic  scattering  cross-section  with  energy  of  electrons 
is  no  longer  negligible  due  to  the  large  number  of 
inelastic  collisions.  To  include  the  energy  dependence  of 
elastic  scattering  cross-section ,  we  use  5L  as  the 
parameter  of  energy  of  an  electron  instead  of  using 
inelastic  scattering  cross-section.  The  energy  of  an 
electron  is  assumed  to  be  a  well  defined  function  of  5L  , 
i.e.  we  assume  "continuous  slowing  down  approximation" 
according  to  the  linear  stopping  power.  In  the  following 
treatment  we  retain  the  same  boundary  condition  as  before. 
To  put  the  equation  (II.  1)  in  a  more  tractable  form,  we 
expand  the  angular  density  function  in  spherical  harmonics; 


oo  m 


f  (*, V,  1)  =  2  ^ 

m=»  n=-m 


(zm+i)  (m-n)/ 
4tt  (rm-rt)l 


Ymn  ( 


(U38) 


and  obtain  from  (II. 1)  and  (11.38) 


H  +  1  Q 


-  n  s  f  xu 

Xp  r 


Ymn  fir) 


x  0“( u-’v',  0  <Jv  <jtr 


(x.3<P 


where 
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cC.  =  S  Y4  5  Yv  Je 


is  a  constant  vector.  The  boundary  conditions  to  be 
satisfied  by  fmn  are 


"fmh  £no  Nfmn(°) 


4  IT 


With  spherical  harmonics  expansion  of  angular  density 
function,  we  easily  recognize  that  the  spatial  distribution 
alone,  regardless  of  the  orientation  of  electrons, 
is  .  The  last  integral  in  (11.39)  can  be  carried 

out  by  expanding  <f  ({Mr,  b  in  Legendre  polynomials,  using 
the  addition  theorem,  orthogonality,  and  normalization  of 
the  spherical  harmonics; 

+  Km-fmn  ~ ^jp  +  ~ p“ ) Tth-m, _  "p')L+i,n-l 
p""j  "Wi,  n-j  ~ (tn+i) (m-v-n— 1)  n— i 

-  (w-*i  +  i)  _  (m+n)  — 1 -o  <X4l) 

7*?  T£ 


where 


Pm(4')|  Axm  T 


(tl-41) 
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Equation  (11.41)  are  a  set  of  coupled  equations  of 
infinite  number.  To  obtain  an  approximate  solution,  we 
must  truncate  ~fvrm  5  at  a  certain  order.  To  estimate  the 
error  involved  in  this  approximation,  we  again  resort  to 
the  distribution  of  the  orientation  of  velocity  regardless 
of  the  spatial  distribution.  By  integrating  (11.41)  over 
all  space,  and  noting  that  the  terms  with  P  0  vanish 
because  of  the  azimuthal  symmetry  of  the  problem,  we  get 


~bVvr\r\  u  tr 
-  +  Krh  rmn 


0 


(XL.  43) 


where  Fmn  iM  ■=■  |  -f^o  ('i?,  X 


The  solution  to  this 


equation  is 


FmdO  = 


2m-hi 

4 


Km  U')  da) 


which  satisfies  the  boundary  conditions,  so  that 


-fdvU)  =  (  -f  C?,v,  1)  dx 

=  ~  2(2W+l)Pm(p)^cp(-(lKwtil)da')  OL.+t) 


which  is  basically  consistent  with  equation  (11.21).  From 
(11.44)  we  find ; 

<(pm(awe))>=  oxp(-j  Kmfs-’)  da')  (l-45) 


For  a  sharply-peaked  distribution  at  small  angle  0  , 


we 
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can  take 


Pm  (v*  e)  * 


-Lm(m+l)  0 

4 


2, 


so  that 


o 


since  the  argument  of  the  exponential  in  equation  (11.45) 
is  relatively  small  for  high  energies(>20  MeV )  for 
small  m  .  Therefore,  to  adequately  approximate  the 
relatively  shallow  depth  region,  we  have  to  truncate  the 
spherical  harmonics  expansion  of  the  angular  distribution 
function,  at  least,  HI  ~  2,/(  .  For  diffuse  angular 
distribution,  equation  (11.44)  obviously  converges  rapidly. 
In  this  case,  truncating  the  spherical  harmonics  expansion 
at  low  would  yield  a  good  approximation.  But,  unlike 
the  case  of  a  peaked  distribution,  it  is  not  easy  to  locate 
a  point  where  this  approximation  becomes  adequate.  Thus, 
we  arbitrarily  take  the  point  at  which  the  average  cosine 
of  the  polar  angle  becomes  €.  '  ,  that  is, 


or 


where  is  the  mass  stopping  power  of 


a  medium. 
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The  average  penetration  depth  of 

calculated  from  (II„46)by 

( 

<2j>  =  1  e  ex1))  <!&' 

■o 


electrons  can  then  be 


We  apply  the  Fourier-Bessel  integral  and  Fourier 
transformation  to  (11.41)  to  find; 


-  3  m-b  n+ 1  -  ( m+-n+ 1)  (m-nrt')  ^ 

■2SL  1  . 

+  (mHiXrm-n-OSnM.rt-ij 

+  iu>  (  (m+Yi+o  +  (i'1''-t->v>3nH,  n  ]  =0  (1.45s) 


with  the  boundary  condition, 


^wniK^>0) 


where 


g»n(k.W.*) 


Oo 

U(at0e;“2Jn(fep)  pdfd? 


In  order  to  investigate  the  electron  transport  in  the 
region  where  the  depth  is  greater  than  the  diffusion  depth 
defined  by  (11.47),  we  retain  only  terms  with 
both  (T\  and  Vd  up  to  (  The  solution  ^oo  Ifc,  hJ  ?  Q,  )  ,  under 

this  approximation  is, 


da.' 


K,U') 


(ujH  bk1) 


0e.5d) 


The  inverse  transformation  of  (11.58)  with  normalization 


gives ; 
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"foe  (P»£>  $0 
where 


1 2-Jr  x5/1  «*P 


24X 


(p^+tz1) 


(l.5l) 


(X.52) 


For  the  numerical  calculation  of  spatial 
distribution  in  the  depth  region,  we  need  the  mass  stopping 
power  and  the  energy  dependence  of  the  scattering 
cross-section .( 43 )  For  simplicity,  the  Rutherford  formula 
accounting  for  inelastic  deflections  is  used.  The 
"screening"  effect  of  atomic  electrons  is  known  to  be 
negligible  for  our  application  (43),  and  is  omitted  in  our 
calculation.  Thus, 

zirn  <r(vvj o  =  ^  (2+i)(Nk<J>.VaKt+i)z 

x  T_i(T+2)‘2  (l-  c©-^)-2  (x.53) 


where 


Nf,  is  the  Avogadro ' s  number , 

A  is  mass  number, 

4>„  =  8TTe4/3nn2c+ 

and  the  functional  relationship  between  1  and  Si  is  given  as 


follows ; 


I 


. 
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0M=(  _ (§cm_i)  (ir.54) 

)t  l-dT/pda| 

where  (“  dT/ P  d  1 )  is  the  mass  stopping  power.  Equation 
(11.54)  hinges  on  the  legitimacy  of  the  "continuous  slowing 
down  approximation" (csda) .  In  this  approximation,  the 
energy  of  an  electron  is  supposed  to  be  determined  solely 
by  the  total  track-length  which  the  electron  has  travelled. 
A  complete  tabulation  of  mass  stopping  power  and  range  data 
can  be  found  in  Berger  and  Seltzer. (44)  Spencer  considers 
the  stopping  power  number  B  which  is  essentially  a 
logarithmic  function  of  the  energy  to  be  a  constant.  In 
this  approximation  the  mass  stopping  power  may  be  written 
as ; 

=  b(i-(t+')'T 

from  which  follows  the  expression  for  (11.62) 

(3Na<{).Z  B/4A)"1  TMT+lf1 


(H..55) 


(E.5b) 


Spencer  also  presented  the  function  form  of  Km  (0.)  . 

K„  m  -  i  I  z  *  0  (N,<t>.*  /  A)  C„  (<211| 

(3T.5b 

where  Cyn’S  satisfy  the  following  recurrence  relations. 


Co  =  o 


75 


c,  = 

Cm+i  ~  (2-+W')Cm  -  O+rfl'1)  Cm-i  "  (24'rn’>) 

If,  therefore,  we  define  c{  -  (^Na^o  2  B/4A)  ^4.  ,  we  may 

rewrite  (11.56)  in  the  form  TY(T+I)  =  I4»-£(t')]-4/o(  . 

| 

Defining  dm  E  ^  (2  +  1 )  Cm/  8  ,  it  is  easily  seen  that 

(11.57)  may  be  represented  with  considerable  accuracy  by 
the  form; 

Kni(JL)sj  ddm/[i.-i(T^]CA.-|L(T)+^] 

The  spatial  distribution  of  electrons  in  the  depth 
region  then  can  be  calculated  from  (11.51). 

(  dv-f  |  f00  (P,  2,  £  )  4S. 

o 


1 1 1 . 4  Discussion 

Assuming  that  $1  represents  the  estimate  of  the 
path-length  of  individual  electrons,  equation  (11.21)  means 
the  orientations  of  the  electrons  are  more  random  for 
larger  or  equivalently,  lower  energy  electrons .( Figure 
14)  In  this  calculation,  as  electrons  are  assumed  to  be 
placed  at  the  surface  of  a  semi-infinite  medium  and 
oriented  perpendicularly  to  the  surface,  it  is  likely  that 
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y  =  cos  0 


Fig.  14.  The  plot  of  the  angular  distribution  function  f (0,t) 
for  various  i  .  f (0,t)  is  normalized  to  the  f (0 ,1) . 

Z  is  the  pathlength  of  an  electron  divided  by  the 
diffusion  length. 


most  of  electrons  with  large  angle  orientation  are  those 
which  have  relatively  large  X  .  The  contribution  from  a 
single  large  angle  scattering  to  the  number  of  electrons  of 
large  angle  orientation  would  likely  be  significant  after 
electrons  has  travelled  for  a  long  time  "X  ">  because  the 
single  large  angle  scattering  process  is  not  a  frequent 
event . 


As  for  Fermi's  model (46)  of  electron  transport, 
since  it  incorporates  the  small  angle  approximation,  it 
underestimates  the  probability  that  an  electron  has  a  large 
angle  orientation ( Figure  15).  This  is  manifested  in  the 
expression  (11.23).  That  is,  the  prediction  obtained  by 
Fermi  equation  differs  from  that  obtained  by  the  Boltzmann 
equation  by  as  much  as  (Al/T\  0/0  )[^z  even  in  the  superficial 
region,  &  .  This  trend  is  expected  to  be  more 
pronounced  at  the  deeper  region,  due  to  the  reason 
explained  previously. 

In  addition  to  the  inadequate  description  of  angular 
distribution,  Fermi's  model  is  not  appropriate  for  the 
description  of  spatial  distribution  of  the  incident 
electrons  in  the  deeper  region  of  semi-infinite  medium 
because  it  does  not  include  any  delimiting  parameter  for 
the  stopping  of  electrons.  This  is  because  this  model  is 
basically  medium-oriented,  that  is,  the  primary  objective 
of  it  is  to  get  the  distribution  of  electrons  after  they 


(e) 1 
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Fig.  15.  The  angular  distribution  of  elections.  (56) 
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have  been  transported  through  a  certain  thickness  of  a 
slab®  The  distribution  of  electrons  after  they  have  lost  a 
certain  amount  of  their  initial  energy,  is  not  considered 
and  is  more  relevant  to  transport  theory.  The  depth  of  a 
medium  is  not  a  suitable  quantity  to  represent  the  physical 
state  of  individual  electron  because  the  fluctuation  of 
energy  loss  and  scattering  spoils  the  relationship  between 
the  longitudinal  position  of  electrons  and  their  real 


physical  state. 


CHAPTER  IV 

CONCLUSION  AND  FUTURE  WORK 


The  Boltzmann  equation  can  be  used  to  obtain  the 
spatial  fluence  distribution  of  electrons  which  has 
immediate  application  to  dose  calculation  in  a  complex 
medium.  In  addition,  this  equation  provides  the 
predictability  of  dose  distribution  without  resorting  to 
experimental  parameters.  Unfortunately,  the  complex  nature 
of  the  Boltzmann  equation  seems  to  defy  a  solution  in 
closed  form,  at  least  for  a  pencil  beam  boundary  condition. 
However,  with  some  reasonable  approximations,  the  equation 
can  be  utilized  to  investigate  some  salient  features  of 
electron  transport  which  would  be  useful  for  seeking 
numerical  solutions,  and  designing  future  experiments. 

Despite  its  potential  impact  on  electron  radiation 
therapy,  the  electron  transport  problem  has  been  largely 
ignored.  Several  misinterpretations  or  incomplete 
interpretations  of  experimental  results  have  originated 
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from  this  ignorance.  For  example,  current  understanding  of 
the  macroscopic  effects  of  inhomogeneities  is  mainly 
related  to  the  density  of  inhomogeneities.  In  fact,  as 
described  in  the  first  chapter,  because  the  microscopic 
elastic  interactions  between  an  incident  electron  and  a 
constituent  atom  is  influenced  by  the  atomic  number,  the 
effect  of  inhomogeneities  should  be  regarded  as  the 
combined  action  of  the  change  of  the  density  and  the  atomic 
number  of  the  medium. 

Regarding  the  atomic  number  effect  of 
inhomogeneities,  it  is  expected  that  the  anisotropy  of 
elastic  scattering  following  the  change  of  the  atomic 
number  would  significantly  contribute  to  the  distortion  of 
electron  fluence  distribution  in  a  homogeneous  medium,  and 
thus  affect  the  dose  distribution.  This  contribution 
should  be  more  significant  while  electrons  are  retaining 
higher  energies,  that  is,  at  shallow  depth  of  the  medium. 
This  is  not  only  because  the  anisotropy  of  elastic 
scattering  cross-section  is  more  pronounced  at  higher 
energies,  but  also  because  at  these  higher  energies,  most 
of  electrons  are  still  well  correlated  spatially.  Thus, 
the  breadown  in  correlation  would  be  more  appreciable. 
This  contribution  should  become  less  pronounced  as  the 
orientations  of  electrons  become  diffuse  at  deeper 
locations,  because  in  this  case,  the  anisotropy  of 
scattering  cross-section  would  not  be  as  influential  as 
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before  in  changing  the  electron  fluence  distribution.  Even 
in  a  medium  which  is  made  up  of  the  atoms  whose  scattering 
cross-sections  are  highly  anisotropic,  a  point  source 
particles  would  ultimately  produce  a  spherically  symmetric 
particle  fluence  distribution. 

Better  understanding  of  the  electron  transport 
problem  can  be  utilized  to  a  great  advantage  for  innovation 
in  the  current  dose  administration  technique.  For  example, 
medical  linear  accelerators  which  use  magnetically  scanned 
electron  beams  have  some  advantageous  features,  described 
in  the  first  chapter,  which  cannot  be  found  in  the 
conventional  accelerators  using  foil-scattering  mechanism. 
However,  to  date,  this  magnet-scanning  technique  has  only 
been  regarded  as  one  of  the  possible  methods  of  producing 
broad  beam  without  perception  of  other  inherent  advantages. 
It  is  conceivable  that  this  technique  would  eliminate,  in 
addition  to  the  scattering  foil,  the  necessity  for; (a) 
field  shaping  devices,  including  the  collimators,  lead 
masks,  and  (b)  external  compensating  materials  for  the 
effect  of  inhomogeneity.  This  technique  can  eliminate  the 
former  by  choosing  an  appropriate  scan  pattern  and  the 
latter  by  modulating  the  intensity  and/or  energy  of  the 
each  narrow  beams.  More  standardization  of  a  practical 
electron  beam  is  also  to  be  expected  by  using  this 


technique . 


' 
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It  is  certain  that  the  essential  ground  work  for  the 
new  generation  of  spot  beam  scanning  technique ( 39 ) ,  should 
be  the  development  of  transport  theory  for  electrons.  Even 
though  the  ultimate  information  we  need  is  the  dose 
distribution,  the  primary  physical  quantity  to  be  sought  in 
radiation  physics  should  be  the  fluence  of  electrons.  In 
this  context,  a  new  device  for  measuring  the  particle 
fluence  with  reasonable  resolution  is  necessary  to 
compliment  the  current  dosimetric  devices,  and  confirm 
theoretical  expectations. 
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APPENDIX  A 

RELATIONSHIP  OF  BOLTZMANN  AND  FERMI  EQUATIONS 

As  a  first  step  in  deriving  the  Fermi  equation  from  the 
linearized  Boltzmann  equation  (II.l),  we  have  to  modify  the 
equation  (II.l)  slightly,  because  there  are  substantial 
differences  between  the  two  models,  notably; 

(a)  In  Fermi's  model,  electrons  remains  on  a  plane 
containing  the  initial  trajectory  after  scattering, 
that  is,  scattering  process  changes  only  the  polar 
angle  of  the  orientation  of  the  incident  electron. 

(b)  As  the  energy  losses  of  electrons  are  ignored  in 
Fermi's  model,  the  path-length  of  an  electron,  ,  is 
not  considered. 

(c)  Fermi's  model  defines  a  probability  that  an  electron 
traversing  the  thickness  At  of  a  scattering  medium 
will  be  deflected  by  an  angle  ( 0  ,  d©  )  regardless  of 
the  detailed  multiple  scattering  processes. 
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Due  to  this  particular  treatment  of  multiple  scattering, 
the  necessary  changes  of  the  equation  (II *1)  for  the 
comparison  are; 

(a)  There  is  no  (fy  -dependence  in  Q~  and  *p  .  In  addition, 
the  scattering  cross-section  becomes  a  function  of 

\\)  =  CfluT1  (u-  ij')  “  O'"  B 

With  these  changes  the  equation  (II. 1)  reads; 

(^0)  -  n(  de'  i -  f  (>?, 0)  o~(©“  ©)  (A  -  0 


(b)  The  ^  -dependence  of  the  angular  distribution 
function,  -f  ,  and  the  scattering  cross-section,  (T  ,  in 
the  equation  (II.l)  should  be  ignored. 

AS  cfe'-e)  is  peaked  at  very  small  (©  — 0"),  we 
expand  "f  (V,  0’)  in  a  Taylor  series  of 

f  (?>  ©')  ~  f  0?,e)  +  ^ 


Thus,  the  right  hand  side  of  (A.l)  becomes; 


n 


[  de'(e'-e) <r(e'-e)j  ~  +  |peTe'-e)We)"§>}_ 


The 


limit  of  integration  can  be 


extended 


,  l 
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from  —  Oo  to  Do  ,  since  (T(B/~0')  is  very  sharply  peaked 
a. t  0  =  0.  Since  O”(0i  -0s)  is  even  in  (O'-Q)  ,  the  equation 
(A.l)  reduces  to; 


+  C&69 


^2 


O© 


n  de'  cr(©'-©")  (e'-eV 


Oo 


Tie* 


According  to  the  Fermi's  treatment  of  angular  deflection, 
explained  in  (c),  the  expression  in  the  bracket  of  the 
equation  (A. 2)  is  equivalent  to  l/tQ1  .  Finally, 
since  0  is  small,  (A. 2)  reduces  to  Fermi  equat ion . ( 49 ) 

-  _0.lf  +  i 

T7z  "^P  Co"  7>d* 


Q .  E  .  D  . 


